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ABSTRACT. We consider the exclusion process in the one-dimensional discrete 
torus with N points, where all the bonds have conductance one, except a finite 
number of slow bonds, with conductance N^ 13 , with £ [0, oo). We prove that 
the time evolution of the empirical density of particles, in the diffusive scaling, 
has a distinct behavior according to the range of the parameter /3. If j3 g [0, 1), 
the hydrodynamic limit is given by the usual heat equation. If /3 = 1, it is 
given by a parabolic equation involving an operator where W is the 

Lebesgue measure on the torus plus the sum of the Dirac measure supported 
on each macroscopic point related to the slow bond. If (3 G (1, oo), it is given by 
the heat equation with Neumann's boundary conditions, meaning no passage 
through the slow bonds in the continuum. 



1. Introduction 

An important subject in statistical physics is the characterization of the 
hydrodynamical behavior of interacting particle systems in random or inho- 
mogeneous media. One relevant and puzzling problem is to consider particle 
systems with slow bonds and to analyze the macroscopic effect on the hydro- 
dynamic profiles, depending on the strength at these bonds. The problem we 
address in this paper is the complete characterization of the hydrodynamic 
limit scenario for the exclusion process with a finite number of slow bonds. De- 
pending on the strength at the slow bonds, one observes a change of behavior 
that goes from smooth profiles to the development of discontinuities. 

We begin by giving a brief and far from complete review about some results 
on the subject, all of them related to the exclusion process. In (3), by taking 
suitable random conductances {ck : k > 1}, such that {c^ 1 : k > 1} satisfy a 
Law of Large Numbers, it was proved that the randomness of the medium does 
not survive in the macroscopic time evolution of the density of particles. In flU, 
the authors consider conductances driven by an a-stable subordinator W, and 
in this case, the randomness survives in the continuum, by replacing in the hy- 
drodynamical equation the usual Laplacian by a generalized operator 
which results in the weak heat equation. In the same line of such quenched 



2000 Mathematics Subject Classification. 60K35,26A24,35K55. 

Key words and phrases. Hydrodynamic limit, exclusion process, slow bonds. 

1 



2 



TERTULIANO FRANCO, PATRfCIA GONCALVES, AND ADRIANA NEUMANN 



result, E5] shows the analogous behavior, but for a general strictly increasing 
function W. All the previous works are restricted to the one-dimensional set- 
ting, and strongly based on convergence results for diffusions or random walks 
in one-dimensional inhomogeneous media, see II 141 . In II15L there is a general- 
ization of |5] for a suitable d-dimensional setting, in some sense decomposable 
into d one-dimensional cases. General sufficient conditions for the hydrody- 
namical limit of exclusion process in inhomogeneous medium were established 
in 0. All the above works have in common the association of the exponen- 
tial clock with the bonds, having the Bernoulli product measure as invariant 
measure, and being close, in some sense, to the symmetric simple exclusion 
process. 

In fl3l . the totally asymmetric simple exclusion process is considered to 
have a single bond with smaller clock parameter. Such "slow bond", not only 
slows down the passage of particles across it, but it also has a macroscopical 
impact since it disturbs the hydrodynamic profile. Somewhat intermediate be- 
tween the symmetric and asymmetric case, in QJ is considered a single asym- 
metric bond in the exclusion process. This unique asymmetric bond gives rise 
to a flux in the torus and also influences the macroscopic evolution of the den- 
sity of particles. In the symmetric case, obtained a d-dimensional result for 
a model in which the slow bonds are close to a smooth surface. 

As a consequence of the above results, one can observe the recurrent phe- 
nomena about the distinct characteristics of slow bonds in symmetric and 
asymmetric settings. In the asymmetric case, e.g. Ifl3ll and (U, the slow bond 
parameter does not need to be rescaled in order to have a macroscopic influ- 
ence. Nevertheless, in the symmetric case, from (4), (BJ and we see that the 
slow bond must have parameter of order iV _1 in order to have macroscopical 
impact. 

In this paper, we make precise this last statement for the following model. 
Consider the state space of configurations with at most one particle per site in 
the discrete torus. To each bond is associated an exponential clock. When this 
clock rings, the occupancies of the sites connected by the bond are exchanged. 
All the bonds have clock parameter equal to 1, except k finite bonds, chosen in 
such a way that these bonds correspond to k fixed macroscopic points b\, . . . , bk. 
The conductances in these slow bonds are given by N~P, with (3 e [0, +oo) and 
the scale here is diffusive in all bonds. 

If j3 = 1, the time evolution of the density of particles p(t, •) is described by 
the partial differential equation 
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where the operator ^ is defined in Subsection 12.11 and W is the Lebesgue 
measure on the torus plus the sum of the Dirac measure in each of the {hi : i = 
1, k}. This result is a particular case of both the results in and 0. For 
the sake of completeness, we present here a simpler proof of it. It is relevant 
to mention the interpretation of such partial differential equation as a weak 
version of 

dtp = dip 
< d uPt (l) = d uPt (0) = pt{l)-p t (0) 
, p(Pr) = 7(0 

where and 1 mean the left and right side of a macroscopic point 6, related to 
a slow bond. This equation says that p is discontinuous at each macroscopic 
point {b t : i = 1, ...,k} with passage of mass at such point and is governed 
by the Fick's Law: the rate of passage of mass is proportional to the gradient 
concentration. Such interpretation comes from the natural domain Cw of the 
operator defined in Subsection 16.21 It is easy to verify that all the func- 

tions in the domain Cw satisfy the above boundary condition, for more details 
see © and 0. 

If P e [0, 1), the conductances in these slow bonds do not converge to zero 
sufficiently fast in order to appear in the hydrodynamical limit. As a conse- 
quence, there is no macroscopical influence of the slow bonds in the continuum 
and we obtain the hydrodynamical equation as the usual heat equation. The 
proof of last result is based on the Replacement Lemma, and the range param- 
eter of P is sharp in the sense that, it only works for f3 e [0, 1). 

As p increases, the conductance at the slow bonds decreases and the pas- 
sage of particles through these bonds becomes more difficult. In fact, for 
P e (1, +00), the clock parameters go to zero faster than at the critical value 
p = 1 and each slow bond gives rise to a barrier in the continuum. Macroscop- 
ically this phenomena gives rise to the usual heat equation with Neumann's 
boundary conditions at each macroscopic point {b % : i = 1, k}, which means 
here that the spatial derivative of p at each {bi : i = 1, fc} equals to zero and, 
physically, this represents an isolated boundary. Moreover, the uniqueness of 
weak solutions of such equation says explicitly that the macroscopic evolution 
of the density of particles is independent for each interval [bi,b i+ {\, however 
the passage of particles in the discrete torus through the slow bonds is still 
possible. The proof of this result is also based on the Replacement Lemma and 
requires sharp energy estimates. 

Since the regime p = 1 was already known from previous works, the main 
contribution of this article is the complete characterization of the three dis- 
tinct behaviors for the time evolution of the empirical density of particles, ex- 
hibiting a behavior change depending on the parameter of the conductance at 
the slow bonds. From our knowledge, no similar phenomena were exploited 
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for the hydrodynamic limit of interacting particle systems. Moreover, for the 
regime /3 e (1, oo) the density evolves according to the heat equation with Neu- 
mann's boundary conditions, which has a meaningful physical interpretation. 
This is the other great novelty developed in this paper. So far, partial differ- 
ential equations with Dirichlet's boundary conditions could be approached by 
e.g. studying interacting particle systems in contact with reservoirs. Here, 
by considering partial differential equations with Neumann's boundary condi- 
tions, we give a step towards extending the set of treatable partial differen- 
tial equations by the hydrodynamic limit theory. Besides all the mentioned 
achievements, we also prove that the regime p = 1 is critical, since the other 
two regimes have positive Lebesgue measure on the line. 

In order to achieve our goal, the main difficulties appear in the character- 
ization of limit points for each regime of /3. We overcome this difficulty by 
developing a suitable Replacement Lemma, which allow us to replace product 
of site occupancies by functions of the empirical measure in the continuum 
limit. Furthermore, that lemma is also crucial for characterizing the behavior 
near the slow bonds. 

Our result can also be extended to non-degenerate exclusion type models as 
introduced in [8j. In such models, particles interact with hard core exclusion 
and the rate of exchange between two consecutive sites is influenced by the 
number of particles in the vicinity of the exchanging sites. The jump rate is 
strictly positive, so that all the configurations are erdogic, in the sense that 
a move to an unoccupied site can always occur. It was shown in O that the 
hydrodynamical equation for such models is given by a non-linear partial equa- 
tion. Having established the Replacement Lemma, the extension of our results 
to these models is almost standard 0. We also believe that our method is 
robust enough fitting other models such as independent random walks, the 
zero-range process, the generalized exclusion process, when a finite number of 
slow bonds is present. 

The present work is divided as follows. In Section [21 we introduce notation 
and state the main result, namely Theorem l2.ll In Section [3] we make precise 
the scaling limit and sketch the proof of Theorem 12. II In Section [H we prove 
tightness for any range of the parameter /?. In Section^ we prove the Replace- 
ment Lemma and we establish the energy estimates, which are fundamental 
for characterizing the limit points and the uniqueness of weak solutions of the 
partial differential equations considered here. In Section[6]we characterize the 
limit points as weak solutions of the corresponding partial differential equa- 
tions. Finally, uniqueness of weak solutions is refereed to Section [71 



EXCLUSION PROCESS WITH SLOW BONDS 



5 



2. Notation and Results 

Let Tjv = {1, ■ ■ • , N} be the one-dimensional discrete torus with ./V points. At 
each site, we allow at most one particle. Therefore, we will be concerned about 
the state space {0, 1} Tjv . Configurations will be denoted by the Greek letter 77, 
so that rj(x) = 1, if the site x is occupied, otherwise rj(x) = 0. 

We define now the exclusion process with state space {0, 1} Tjv and with con- 
ductance {^x,x+x}x at the bond of vertices x,x+l. The dynamics of this Markov 
process can be described as follows. To each bond of vertices x, x + 1, we as- 
sociate an exponential clock of parameter C^+i- When this clock rings, the 
value of rj at the vertices of this bond are exchanged. This process can also be 
characterized in terms of its infinitesimal generator Cn, which acts on local 
functions / : {0, 1} T " -> K as 

where rf' x+1 is the configuration obtained from r\ by exchanging the variables 

T](x) and t](x + 1): 

77(21 + 1), if y = x , 
'■ X+1 )(V)={ ri{x), ify = x + l, 
r](y), otherwise. 



The Bernoulli product measures : < a < 1} are invariant and in 
fact, reversible, for the dynamics introduced above. Namely, u£ is a product 
measure on {0, 1} T ™ with marginal at site x in Tn given by 

v a {n : V(x) = 1} = a. 

Denote by T the one-dimensional continuous torus [0, 1). The exclusion pro- 
cess with a slow bond at each point b± . . . , b k e T is defined with the following 
conductances: 



N-f>, if {h,...,b k }n(§ ,^±1] ^0, 



1, otherwise. 



The conductances are chosen in such a way that particles cross bonds at 
rate one, except k particular bonds in which the dynamics is slowed down by 
a factor N~ fj , with @ e [0, 00). Each one of these particular bonds contains the 
macroscopic point 64 e T; or 6, coincides with some vertex and the slow bond 
is chosen as the bond to the left of ^ . To simplify notation, we denote by Nbi 
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the left vertex of the slow bond containing b t . 

Denote by {?/ t := rj tN 2 : t > 0} the Markov process on {0, 1} Tn associated to 
the generator £ N speeded up by TV 2 . Although rjt depends on N and j3, we are 
not indexing it on that in order not to overload notation. Let D(R + , {0, 1} Tjv ) be 
the path space of cadlag trajectories with values in {0, l} Tlv . For a measure /z/v 
on {0, l} TiV , denote by P^ v the probability measure on D(R+, {0, 1} Tjv ) induced 
by the initial state /jn and the Markov process {rj t : t > 0} and denote by E^ N 
the expectation with respect to P^ N . 



Definition 1. A sequence of probability measures {^n ■ N > 1} on {0, 1} Tjv 
is said to be associated to a profile po ■ T — > [0,1] if for every 8 > and every 
continuous functions H : T — > K 

lim m U: U H(§) V (x)- f H(u)p (u)du >S\ = 0. (1) 

Now we introduce an operator which corresponds to the generator of the 
random walk in T N with conductance ^x+i a ^ the bond of vertices x, x + 1. 
This operator acts on H : T — >• R as 



(2) 



We will not differentiate the notation for functions H denned on T and on 
T N . The indicator function of a set A will be written by 1a{u), which is one 
when u £ A and zero otherwise. 

2.1. The Operator A 

Given the points b\, . . . , bk G T, define the measure W(du) in the torus T by 

W(du) = du + S bl (du) + ■ ■ ■ + Sb k (du) , 

so that W is the Lebesgue measure on the torus T plus the sum of the Dirac 
measure in each of the {&; :i = l, k}. 

Let Hly be the set of functions F in L 2 (T) such that for x £ T 



F(x) = a + 6+ / f(z)dz)W(dy), 

J(0,x] V JO 

for some function / in L 2 (T) and a, b e R such that 



/(.T)d.T = 0, / (b+ f{z)dz)w(dy) = 
J(o,i] Jo 7 



(3) 
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Define the operator 

JL—f = f. 

dx dW J 
For more details we refer the reader to |5|. 

2.2. The hydrodynamical equations. 

Consider a continuous density profile 7 : T — > [0, 1]. Denote by (•, •) the inner 
product in L 2 (Y), by p t a function p(t, •) and for an integer n denote by C n (T) 
the set of continuous functions from T to R and with continuous derivatives of 
order up to n. For 1 an interval of T, here and in the sequel, for n and m inte- 
gers, we use the notation C™' m ([0, T] x 1) to denote the set of functions defined 
on the domain [0, T] x I, that are of class C" in time and C m in space. 



Definition 2. A bounded function p : [0, T] x T — > K is said to be a weak solution 
of the parabolic differential equation with initial condition j(-): 

S d *P = d uP (4) 

\p(O,0=7(0 

if, for t e [0,T] and H e C 2 (T), p(t, •) satisfies the integral equation 

(pt,H) - (>y,H) - / (Ps,d 2 u H)ds =0. 
Jo 

Definition 3. A bounded function p : [0, T] x T -> K is said to oe a w;ea£ solution 
of the parabolic differential equation with initial condition ■y(-): 

» d d 

dtP = TxdW P (5) 
. P(O,0 = 7(0 

i/J /or £ e [0, T] and H e "H^, •) satisfies the integral equation 

( ft ,a)-(7,a>-jf( A .|^H>4-o. 

Following the notation of [2], denote by L 2 (0, T; 'H 1 (a, 6)) the space of func- 
tions q e £ 2 ([0, T] x [a, b]) for which there exists a function in L 2 ([0, T] x [a, 6]), 
denoted by d u g, satisfying 

nT rb j-T rb 

I I (d u H)(s,u) g{s 7 u) du ds = — / H(s,u) (d u g)(s,u) du ds , 

JO J a JO Ja 

for any H e C 01 ([0, T] x [a, 6]) with compact support in [0, T] x (a, 6). 
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Definition 4. Let [6», 6»+i] C T. A bounded function p : [0,T] x [b» , 6»+i] -> K is 
said to oe a wea& solution of the parabolic differential equation with Neumann's 
boundary conditions in the cylinder [0, T] x [6f, 6i+i] and with initial condition 

7(0: 

ftp = 

p(O,-)=7(0 (6) 
[ d u p(t,bi) = d uP (t,b i+1 ) = 0, Vi G [0,T] 

i/J /or f G [0, T] and if € C 1,2 ([0, T] x [6,, •) satisfies the integral equa- 

tion 

bi+i fbi+i 

p(t,u) H(t,u) du — I 7(w) H(0,u) du 

bi Jbi 
t rbi+i 

p(s, u) {dlH(s, u) + d s H(s, u)} du ds (7) 



Jb, 

+ f d u H{s,b. l+1 )p{ S ,b- +l )ds- f d u H{s,bi)p{sM) d s = 
Jo Jo 

and p(t, •) belongs to L 2 (0, T; b l+1 )). 

Since in Definition |4] we impose p G L 2 (Q, T; % 1 (b l , 6 i+1 )), the integrals are 
well-defined at the boundary. This is a consequence of the following two facts. 
On one hand, it follows from the assumption that p(t, •) e V}{bi, b l+ i), almost 
surely in t G [0, T]. On the other hand, it is well-known that functions belong- 
ing to Ti l (bi,bi + i) and with sided limits at h and b i+ \ are absolutely continu- 
ous with respect to the Lebesgue measure, see Ifl2l for instance. We refer the 
reader to [2] for classical results about Sobolev spaces. 

Heuristically, in order to establish an integral equation for the weak solu- 
tion of the heat equation with Neumann's boundary conditions as above, one 
should multiply I© by a test function H and perform twice a formal integration 
by parts to arrive at ©. 

We are now in position to state the main result of this paper: 

Theorem 2.1. Fix (3 G [0, oo). Consider the exclusion process with k slow bonds 
corresponding to macroscopic points bi,...,bk £ T and with conductance N^ 13 
at each one of these slow bonds. 

Fix a continuous initial profile 7 : T — > [0, 1]. Let {/xjv . N > 1} be a sequence 
of probability measures on {0, 1} Tn associated to 7. Then, for any t G [0, T], for 
every 5 > and every H G C(T), it holds that 

to^K^I- : \t? E H(§)th(x)- J H(u)p(t,u)du\ >S) = 0, 

where : 
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• if P G [0, 1), p(t, ■) is the unique weak solution of ©; 

• if (3 = I, p(t, •) is the unique weak solution of ©; 

• if P £ (1, oo ), in each cylinder [0, T] x [b i: b i+1 ], p(t, •) is the unique weak 
solution of ©. 

Remark 2.2. The assumption that all slow bonds have exactly the same con- 
ductance is not necessary at all. In fact, last result is true when considering 
each slow bond containing the macroscopic point bi with conductance N~°*. In 
that case, we would obtain a parabolic differential equation with the behavior 
at each [bi, &i+i] given by the regime of the corresponding Pi as above. Another 
straightforward generalization is to consider conductances not exactly equal to 
N~P, but of order N" 13 , in the sense that the quotient with N~ f3 converges to 
one. For sake of clarity, we present the proof under the conditions of Theorem 

3. Scaling Limit 

Let M be the space of positive measures on T with total mass bounded by 
one, endowed with the weak topology. Let £ M be the empirical measure 
at time t associated to r\ t , namely, it is the measure on T obtained by rescaling 
space by N and by assigning mass A^ 1 to each particle: 

*f = T E nt{x)S x/N , (8) 

where S u is the Dirac measure concentrated on u. For an integrable function 
H : T -> E, (71-f , H) stands for the integral of H with respect to vrf : 

This notation is not to be mistaken with the inner product in L 2 (R). Also, when 
Ti t has a density p, namely when 7r(£, du) = p(t, u)du, we sometimes write (pt,H) 
for (n t ,H). 

Fix T > 0. Let D([0,T],M) be the space of .M-valued cadlag trajectories n : 
[0, T] -> M endowed with the Skorohod topology For each probability measure 
UN on {0,1} T ™, denote by the measure on the path space D([0,T],M) 
induced by the measure pn and the empirical process ir^ introduced in ®. 

Fix a continuous profile 7:!-) [0,1] and consider a sequence {pn ■ N > 1} 
of measures on {0, 1} Tjv associated to 7. Let Q' 3 be the probability measure on 
D([0, T],M) concentrated on the deterministic path n(t, du) = p(t, u)du, where: 

• if P G [0, 1), p(t, •) is the unique weak solution of ©; 

• if P = 1, p(t, •) is the unique weak solution of ©; 
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• if /3 G (l,oo), in each cylinder [0, T] x [bi, 6»+i], p(t, •) is the unique weak 
solution of©. 

Proposition 3.1. As N t 00, £/ie sequence of probability measures ■ N > 

1} converges weakly to Q 13 . 

The proof of this result is divided into three parts. In the next section, we 
show that the sequence {Q^ : N > 1} is tight, for any j3 G [0, 00). In Section[6] 
we characterize the limit points of this sequence for each regime of the param- 
eter j3. Uniqueness of weak solutions is presented in Section[7]and this implies 
the uniqueness of limit points of the sequence ■ N > 1}. In the fifth 

section, we prove a suitable Replacement Lemma for each regime of /3, which 
is crucial in the task of characterizing limit points and uniqueness. 

4. Tightness 

Proposition 4.1. For any fixed f3 G [0, 00), the sequence of measures : 
N > 1} is tight in the Skorohod topology ofD([0, T],M). 

Proof. In order to prove tightness of {71-^ : < t < T} it is enough to show 
tightness of the real-valued processes {(nf,H) : < t < T} for H G C(T). In 
fact, c.f. IfTOl it is enough to show tightness of { (tt^ , H) : < t < T} for a dense 
set of functions in C(T) with respect to the uniform topology. For that purpose, 
fix H G C 2 (T). By Dynkin's formula, 

M t N (H ) = (vrf , I?) - (7^, If) - / iV 2 £jv(7rf , if) <fa , (9) 

Jo 

is a martingale with respect to the natural filtration Tt := er(?7 s : s < t). In or- 
der to prove tightness of {{tt^, H) : N > 1}, we prove tightness of the sequence 
of the martingales and the integral terms in the decomposition above. We start 
by the former. 

We begin by showing that the L 2 (P^ JV )-norm of the martingale above van- 
ishes as TV -> +00. The quadratic variation of (H) is given by 



(M N (H))t = I £ (Vs(x)-Vs(x + l))(H(^)-H(§)) 



2 



ds. (10) 



It is easy to show that (M N (H)) t < j^||d u lf H^,. Here and in the sequel we use 
the notation ||if||oo : = su P«eT \H{u)\. 

Thus, M^{H) converges to zero as N — >• +00 in L 2 (P^ N ). Notice that above 
we used the trivial bound C^+i < 1- By Doob's inequality, for every 5 > 0, 

= 0, (11) 



lim P£ N 



sup \Mt\H)\ > 6 

0<t<T 



llll IT,, SUP liVj 

JV 

which implies tightness of the sequence of martingales {M^{H); N > 1}. Now, 
we need to examine tightness of the integral term in ©. 



EXCLUSION PROCESS WITH SLOW BONDS 



11 



Denote by Tat the subset of sites i e T« such that x has some adjacent 
slow bond, namely, ^ x+1 = N~ fj or x = N'?. The term N 2 C N {^ ' ,H) 
appearing inside the time integral in © is explicitly given by 



N J2 V s (x){H(^) + H(^)-2H(§) 

x£F N 



By Taylor expansion on H, the absolute value of the first sum above is bounded 
by ||c^i/||oo. Since there are at most 2k elements in T N , £ XjX+ i < 1 and since 
there is only one particle per site, the absolute value of the second sum above 
is bounded by 2 fc||9 u i?'|| 0o . Therefore, there exists a constant C := C'(H, k) > 0, 
such that \N 2 £ n (tt^ ,H)\ < C, which yields 



< C\t- 



By Proposition 4.1.6 of iflOl . last inequality implies tightness of the integral 
term. This concludes the proof. □ 



5. Replacement Lemma and Energy Estimates 

In this section, we obtain fundamental results that allow us to replace the 
mean occupation of a site by the mean density of particles in a small macro- 
scopic box around this site. This result implies that the limit trajectories must 
belong to some Sobolev space, this will be clear later. Before proceeding we 
introduce some tools that we use in the sequel. 

Denote by H N (^, N \v a ) the entropy of a probability measure fi N with respect 
to the invariant state v a . For a precise definition and properties of the entropy, 
we refer the reader to 1101 . In Proposition 18.11 in the Appendix we review a 
classical result saying that there exists a finite constant Kq := Ko(a), such 
that 

H N (n N \v a ) < K Q N, (12) 

for any probability measure fiN £ {0, 1} T,V . 

Denote by (•, -) Va the scalar product of L 2 (v a ) and denote by 33 at the Dirichlet 
form, which is the convex and lower semicontinuous functional (see Corollary 
Al.10.3 of OH) defined as: 



2M/) = (-l n ^1, v7> 
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where / is a probability density with respect to v a (i.e. / > and / fdv a = 1). 
An elementary computation shows that 

fN r 2 

x£T N 

By Theorem Al.9.2 of IflOll . if {Sf : i > 0} stands for the semi-group associated 
to the generator N 2 £ N , then 



H N (jjL N S?\v a ) + N z / Q N (f?)ds < H N (fx N \p a ), 
Jo 

provided stands for the Radon-Nikodym derivative of unS^ (the distribu- 
tion of f] s starting from (U y) with respect to v a . 



5.1. Replacement Lemma. 

Now, we define the local density of particles, which corresponds to the mean 
occupation in a box around a given site. We represent this empirical density in 
the box of size £ around a given site x by j/(x). For e [0, 1), this box can be 
chosen in the usual way, but for j3 e [1, oo), this box must avoid the slow bond. 
From this point on, we denote the integer part of eN, namely L^^J , simply by 
sN. 

Definition 5. For (3 e [0,1), define the empirical density by 

x+eN 
y=x+l 

Definition 6. For (3 G [1, oo), ifx is such that {Nbi, . . . , Nb k }n{x, . . . , x+eN} = 
0, then the empirical density is defined by 

x+eN 

V eN M = E rtv)- 

y—x-\-l 

Otherwise, if, let us say, Nbi G {x, . . . ,x + eN} for some i = 1, .., k, then the 
empirical density is defined by 

Nbi 

V eN (x) = iiv E 

y=Nbi-eN+l 

Since we are considering a finite number of slow bonds, the distance between 
two consecutive macroscopic points related to two consecutive slow bonds is at 
least e, for e sufficiently small. As a consequence, we can suppose, without lost 
of generality that in the previous definition, bi is unique. 
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Lemma 5.1. Fix /3 e [0, 1). Let f be a density with respect to the invariant 
measure v a . Then, 

\v(*)-rf N (x)}f(ri)v a (dTi) < 2(kNP- 1 +e) + NV N (f) 1 VxeT N . 
Proof. From Definition [5] we have that 

x+eN 

[{r,(x)-r 1 sN (x)}f(r ] )u a (d V ) = J £ (v(^) ' v(y))}.f(v) ^(d V ) ■ 
■* ■* y=x+l 

Writing n(x) - n(y) as a telescopic sum, the last expression becomes equal to 

. x+eN y-l 

/{ilv E 52(v(z)-V(z + l))}m" a (dr))- 

y—x+l z—x 

Rewriting the expression above as twice the half and making the transforma- 
tion rj i-> r/ z ' z+1 (for which the probability v a is invariant) it becomes as: 

x+eN y-l 

2m E {r)(z)-v(z+l)}(f(ri)-f(v z ' z+1 ))v a (dv)- 

y—x+l z—x 

Since (a - b) = (y/a - \fb){y/a + Vb) and by the Cauchy-Schwarz's inequality, 
for any A > 0, we bound the previous expression from above by 

x+eN y-l A r , x 9 

E T.TN— / Mz) - V (z + l)} 2 + Vf(v z ' z+1 )) " a (dv) 

y=x+l z=x 'sz.z+l ■> 
x+eN y-l 



t . F JV r / \ 2 

h E E^f 1 J (Vm-VfUP^ Mdn). 



2eN 

y—x+l z—x 

The second sum above is bounded by 

x+eN p 2 

2iN E E ^f- / [Vm- s/JW^j) ".w = &w<j). 

y=a+l 2 eTjv 

On the other hand, since / is a density, the first sum is bounded from above by 

x+eN y-l . . x+eN 

M E Et^^ E 2A(kN^+eN)=2A(kN^+eN). 

y = x + l z = x ^z,z + l y=x + l 

Notice that the term kN^ comes from the existence of k slow bonds. Choosing 
A = jj, the proof ends. □ 

Lemma 5.2 (Replacement Lemma). 

Fix (3 e [0, 1). Let b e T and let x be the right (or left) vertex of the bond 
containing the macroscopic point b. Then, 



:logE„„ 



exp{ 7 iv|^ {r, s (x)-r,f( X )}ds\}] 



(13) 
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Proof. From Jensen's inequality together with the entropy inequality (see for 
example Appendix 1 of iflOll ). for any 7 e R (which will be chosen large), the 
expectation appearing on the statement of the Lemma is bounded from above 
by 

H N (lJ, N \l/ a ) _\_ 

7 TV jN 

By Proposition 18.11 H N ((i N \is a ) < K N, so that it remains to focus on the 
second summand above. Since < e x + e~ x and 

lim i log(a7v + b N ) = max j Urn ± log a N , lim ± log 6jy} , (14) 

we can remove the modulus inside the exponential. By Feynman-Kac's for- 
mula, see Lemma Al.7.2 of II 101 and Proposition I8.2L the second term on the 
right hand side of ( TT3l is less than or equal to 

t sup { [Mx)-rf"(x)}f(v> a (dTi)-NZ) N (f)}. 

f density ^ J J 

Applying Lemma lSTTl and recalling that 7 is arbitrarily large, the proof finishes. 

□ 

The next two results are concerned with both cases /3 = 1 and /3 e (1, 00). 

Lemma 5.3. Fix (3 e [1, 00). Let f be a density with respect to the invariant 
measure v a . Then, 

{r)(x) - r, sN (x)}f( V )u a (dr]) < N® N {f) + 4e ,Vx e T N . 

Moreover, given a function H : T — > K: 

7t E / H(f){r,(x)- V sN (x)}f(r 1 )u a (dr 1 )<ND N (f) + ^ £ (#(#)) 2 . 

Proof. Recall the Definition [6] Let first x be a site such that there is no slow 
bond connecting two sites in {x, . . . , x + eN}. In this case, 



x-f-eiv 



x+eN 
V=x+1 

and following the same arguments as in Lemma 15.11 we bound the previous 
expression from above by 

x „ x+eN y-l . . 
{H{% )f / A 



EE/ jw-Mz) v(* + mHy/m + v / /(^^ +i )) 2 ^(^) 

— r._L1 , — 5z.2 + l 



2sV 

y—x+1 z—x ^z.z-\-l 



2 ^ E E y ^w*) - ^ + i)} a (v?w - vTm) ^(^), 
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Since i^,z+i = 1 f° r a\\ z e {x, . . . ,x + eN — 1}, it yields the boundedness of the 
previous expression by 

2 ®N(f) 



2eNA 



(#(#))' 



A 

Let now a; be a site such that Nb t e {x,. . . ,x + eN} for some i = 1, . . . , k. In 
this case, 

J H{^){n{x)-t N {x)}f{ri)v a {dri) 

. Nbi (15) 

J y=Nbi-eN+l 

Now we split the last summation into two cases, y > x and y < x and then we 
proceed by writing rj(x) — rj(y) as a telescopic sum as in Lemma IBTT1 Then, by 
the same arguments of Lemma 15. II and since £^ 2 +i = 1 for all z in the range 
{Nb l — eN + 1, . . . , Nb t — 1}, we bound the previous expression by 

2 3M/) 



AeNA(H(§ 



A 



Now the first claim of the lemma follows by taking the particular case H (-^) = 
1 and choosing A= jj. 

Finally, if in dl5D we sum over x e Tjy and then divide by N, one concludes 
the second claim of the lemma. □ 

Lemma 5.4 (Replacement Lemma). 

Fix P G [1, oo). Then, for every x eT N 

Hm lim E^ [ f\ Vs (x) - r, E s N (x)} ds 

e— s-0 N— >oo ^ L Jq 

Moreover, given a function H : T — s- K satisfying 
aZso holds 



IV->oo 



lim lim E 

£->0 N->-oo 



t 





[ w E #(#){%(*) -ryf^)}^ 

JO IT 



0. 



Proof. The proof follows exactly the same arguments in Lemma|5Jl Therefore, 
is sufficient to show that the expressions 

t sup { [{ V (x)-r, eN (x)}f(r))dv a -NV N (f)} 

f density ^ J * 

and 

t SU P { / ]?J2 H (&Wx)-V sN (x)}mdv a -NV N (f)}, 

j density ^ J 

vanish as N — > +00, which is an immediate consequence of Lemma l5.31 □ 
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In the next subsection, we will need the following variation of Lemma l5.3l 

Lemma 5.5. Let H : T — > R and let f be a density with respect to v a . Then, for 
every igTjv 



^ E H(§){v(x)-r,(x + eN)}f( V )v a (dr,) 
<NV N (f) + ^ E (H(f)) 2 {e + N^J2l[ b ,M + e)(f)} 



The proof of the last lemma follows the same steps as above and for that 
reason will be omitted. Nevertheless, we sketch the idea of the proof. One 
begins by writing - rj{x + sN) as a telescopic sum and proceeding as in 
Lemma [5.31 The only relevant difference in this case is that is not possible to 
avoid the slow bonds inside the telescopic sum, and therefore the upper bound 
depends on /3. 

5.2. Energy Estimates. 

We prove in this subsection that any limit point Q* of the sequence {Qf ^ : 
N > 1} is concentrated on trajectories p(t,u)du with finite energy, meaning 
that p(t,u) belongs to some Sobolev space. For (3 e [0, 1), this result is an im- 
mediate consequence of the uniqueness of weak solutions of the heat equation. 
The case j3 = 1 is a particular case of the one considered in (5). Therefore, we 
will treat here the remaining case j3 € (1, oo). Such result will play an impor- 
tant role in the uniqueness of weak solutions of©. 

Let Q* be a limit point of {Q£^ : N > 1} and assume without lost of gener- 
ality that the whole sequence converges weakly to Qj. 

Proposition 5.6. The measure is concentrated on paths n(t, u) = p(t, u)du. 
Moreover, there exists a function in L 2 ([0, T] x T), denoted by d u p, such that 

I (d u H)(s,u) p(s,u) du ds = — / H(s,u) (d u p)(s,u) du ds , 
o Jt Jo Jt 

for all H in C o,1 ([0, T] x T) whose support is contained in [0, T] x (T\{6i, . . . , &&}). 

The previous result follows from the next lemma. Recall the definition of 
the constant Kq given in ( fT2l . 

Lemma 5.7. 



E, 



sup < / / (d u H)(s,u) p(s,u) du ds 
h Wo Jt 

pT ? 



2 



/o Jt 



< K 
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where the supremum is carried over all functions H in C°' 1 ([0, T] x T) with 
support contained in [0, T] x (T\{6i, . . . , bk})- 

We start by showing Proposition 15.61 assuming the last result. Later and 
independently we will prove the previous lemma. 

Proof of Proposition^^ Denote by £ : C ' 1 ( [0, T] x T) -> K the linear functional 
defined by 

£{H) = I [ (d u H)(s,u) p(s,u) duds. 
Jo Jt 

Since the set of functions H e C o,1 ([0, T] x T) with support contained in [0, T] x 
(T\{6i, . . . , b k }) is dense in L 2 ([0, T] x T) and since by Lemmal5771 1 is a Qf-a.s. 
bounded functional in C 0,1 ([0, T] x T), we can extend it to a Q*-a.s. bounded 
functional in L 2 ([0, T] x T). In particular, by the Riesz Representation Theorem, 
there exists a function G in L 2 ([0, T] x T) such that 

£(H) = - [ ( H(s,u)G(s,u) duds . 



This finishes the proof. □ 

For a smooth function H: T — s- R, s > and a positive integer AT, define 

V N (e,H, V ) by 

In order to prove the Lemma [5.7[ we need the following technical result: 

Lemma 5.8. Consider Hi, . . . , H k functions in C 01 ([0, T] x T) with support con- 
tained in [0, T] x (T\{6i, . . . , bk}). Hence, for every e > 0: 

< K . (16) 



lim lim E* 3 „ 

{^0 «^oo f 



max 

Ki<k 



V N (e,H t (8,-),rf a N )d8} 



Proof. It follows from Lemma [5~4l that in order to prove ( [TBI , we just need to 
show that 



lim E P 



max 

Ki<k 



V N {e,Hi(s,-),r] s )ds} 



o ■ 



— log E fc 



exp ■ 



max TV / VKr(e, J?»(s, •), %)ds [■ 

l KKfc /n - 1 



AT-»-oo ^ 

By the entropy and the Jensen's inequality, for each fixed N, the previous ex- 
pectation is less than or equal to 

TV N 

By ( TT2l ). the first term above is bounded by K . Since exp{maxi<j<fe a,} is 
bounded from above by J2i< 3 <k exp{a,} and by ( fl4l ). the limit as TV | oo, of the 
second term of the previous expression is less than or equal to 

max lim -^-logE^ exp{/V~ / Vn{s, H(s, ■), r] s )ds\ 

Ki<k N—t-oo I\ L I In. ) 
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We now prove that, for each fixed i the limit above is nonpositive. 

Fix 1 < i < k. By the Feynman-Kac's formula and the variational formula 
for the largest eigenvalue of a symmetric operator, for each fixed N, the previ- 
ous expectation is bounded from above by 



sup 

/ 



V N (e,Hi(s, •), v)f(v)Mdv) - NV N (f)} 



ds. 



In last formula the supremum is taken over all probability densities / with 
respect to v a . By assumption, each of the functions {Hi : i = 1, . . . , k} vanishes 
in a neighborhood of each 6; e T. This together with Lemma l531 imply that the 
previous expression has nonpositive limsup. This is enough to conclude. □ 

We define now an approximation of the identity in the continuous torus 
given by 

i l {v , v+£) (u) , if v € T\ U.ti (h - e, h) , 
i £ (u,v)=i \ l(ii- £ i,)(i) , if v e (6i - e, bi) , 

i l(b k -e,b k )(u) , if v e {b k - e, b k ) . 
The convolution of a measure tt with l £ is defined by 

(ir*L e )(v) = / L e (u, v) -K(du) . 



(17) 



For a function p, the convolution p * l £ is understood as the convolution of 
the measure p{u)du with l e . Recall Definition [6l At this point, an important 
remark is the equality 

r ) f(x) = (^*L E )(§), (18) 
which is of straightforward verification. 

Proof of Lemma \577\ Consider a sequence {H l : i > 1} dense (with respect to 
the norm + ||<9„if ||oo) in the subset of ^([O.T] x T) of functions with 

support contained in [0, T] x (T\{&i, . . . , bk}). 



Recall that we suppose that 
, for every k > 1, 



,N 



N > 1} converges to Q%. By <[16]> and 



lim E n n 



max (i f [ H % {s,u) { p s s (u) - p s s {u + e)\ du , 

Kl<k I £ l J T 



JT 



(Hi(s, u)) 2 du iis| 



< K , 
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where p 5 s (u) = (p s * ts)(u) as denned above. Letting S 10, performing a change 
of variables and then letting s I 0, we obtain that 

■En, 

/T 



max / / (d u Hi)(s,u)p(s,u) duds 



- 2 / / u)r duds 



< if n 



To conclude the proof it remains to apply the Monotone Convergence Theorem 
and recall that {H l : i > 1} is a dense sequence (with respect to the norm 
Halloo + |j<9 u iJ||oo) in the subset of functions of C o,1 ([0,T] x T) with support 
contained in [0, T] x (T\{6 X . . . , b k }). □ 

Remark 5.9. In terms of Sobolev spaces, we have just proved that, for (3 e 
(1, oo), Q^-almost surely, the limit trajectory p(t, u)du is such that p(t, u) belongs 
to L X (Q, T; H 1 ^, b i+ i)), in each cylinder [0, T] x (b t ,b i+ i). Notice that in view of 
the presence of slow bonds and of Lemma 15.51 is it not possible to obtain the 
same result considering the whole space i 1 (0, T; 'H 1 (T)). 

6. Characterization of Limit Points 

We prove in this section that all limit points Q* of the sequence {Qf ^ : N > 
1} are concentrated on trajectories of measures absolutely continuous with re- 
spect to the Lebesgue measure: n(t,du) = p(t,u)du, whose density p(t,u) is a 
weak solution of the hydrodynamic equation ©, © or ©, for each correspond- 
ing value of p. 

Let Q* be a limit point of the sequence {Qf^ : N > 1} and assume, without 
lost of generality, that : > 1} converges to Q;. The existence of Q; is 

guaranteed by Proposition [4J] 

Since there is at most one particle per site, it is easy to show that Q* is 
concentrated on trajectories 7r t (du) which are absolutely continuous with re- 
spect to the Lebesgue measure, ir t {du) = p(t,u)du and whose density p(-)t, ■ is 
non-negative and bounded by 1 (for more details see IflOl ). We distinguish the 
regime of /3 in different subsections below. In all the cases, we will make use 
of the martingale Mf 1 '(H) defined in I©. By a simple change of variables, the 
integral term in © can be rewritten as a function of the empirical measure, 
such that: 

M t N (H) = (Trf.fl) - (ir^H) - / (vrf , N 2 h N H) ds , (19) 

Jo 

where L w was defined in ©. 

We notice here that, for any choice of H, M^(H) is a martingale. In due 
course we impose extra conditions on H in order to identify the density p(t, ■) as 
a weak solution of the corresponding weak equation depending on the regime 
of the parameter p. 
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6.1. Characterization of Limit Points for [3 e [0, 1). 

Here, we want to show that p(t, ■) is a weak solution of ©. Let H e C 2 (T). 
We begin by claiming that 



7T. : (ir t ,H) - (n ,H) - / {-K s ,d 2 u H)ds = 0, Vt e [0,T] 



= 1. 



(20) 



In order to prove the last claim, it is enough to show that, for every 5 > 0: 



7T. : sup 

0<i<T 



[n t ,H) - (n Q ,H) 



(ir s ,dlH) ds 



> 5 



0. 



By Portmanteau's Theorem and Proposition 18.31 last probability is bounded 
from above by 



lim 0^5 



7T. : sup 

0<t<T 



[w u H) - ( no ,H) - I (n a ,%H)d3 
'o 



> 5 



since the supremum above is a continuous function in the Skorohod metric. 
Adding and subtracting (tt^ 1 , N 2 h^H) in the integral term above and recalling 
the definition of the previous expression is bounded from above by 



lim 



lim 

N— s-oc 



sup 

0<t<T 



sup 

0<i<T 



(< , N 2 h N H) ds > 6/2 



?rf , d 2 u H - N 2 L N H) ds > 5/2 



By dl9l > and 01H . the first term in last expression is null. By the definition of 
T N given in Section |4] and since there is only one particle per site, the second 
term in last expression becomes bounded by 



lim 

N— >oo 



+ lim 

AT— >oc 



N 



> 5/4 



SU P I / W E {d 2 u H(f)-N 2 h N H(f)}r, s (x)ds > 5/4 
o<t<T I Jo x ^ n I- 



Outside r^v, the operator N 2 h N coincides with the discrete Laplacian and 
since H e C 2 (T), the first term in last expression is zero. Recall that there 
are 2k elements in Tn- Applying the triangular inequality, the second expres- 
sion in the previous sum becomes bounded by 



lim 



N- 



lim 



^fiKmu > s/s 



sup I Yl / NL nH(§ )r )s (x)ds > 5/8 



o<t<T ' — Jo 



For large N, the first probability vanishes. Now we deal with the second term. 
We associate to each slow bond containing a point a unique pair of sites in 
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Tn, namely Nbi and Nbi + 1. By the triangular inequality, in order to show 
that the second expression above is zero, it is sufficient to verify that 



lim Pf? 



sup 


L 


L 0<t<T 





{NL N H(Z&)T, a (Nbi) 



-Nh N H(^^-)?] s (Nbi + l)}ds > S/8k 



0. 



l,...,k. The expression inside the integral above can be explicitly 



for each i = 
written as 

{N [if(^pi) - H(^)) + [H(^±l) H^)]}r, a (Nbi) 
+ {N^ [H(^) - H(^tl)} + N [tf (M) _ H(^tl)} } n.iNbi + 1) . 

Since H is smooth and /3 e [0, 1), the terms inside the parenthesis involving 
iV 1- ' 3 converge to zero and the terms involving N converge to plus or minus 
the space derivative of H at hi. Therefore, again by the triangular inequality, 
it remains to show that, for any 5 > 0, 



lim 

JV->oc 



J/3 



sup 


I 


L 0<t<T 





ds 



> 5 



(21) 



equals to zero. The integral inside the probability above is continuous as 
a function of the time t. Moreover, it has a Lipschitz constant bounded by 
\d u H(bi)\. If d u H(bi) = 0, then there is nothing to do. Otherwise, let t = < 
ti < ■ ■ ■ < t n = T be a partition of [0, T] with mesh bounded by 6{\2d u H 
Notice the partition is fixed, depending only on the function H. By the trian- 
gular inequality, d2ll ) is bounded by 



Elim 
AT->o 



3=0 



d u H(bi){ri s (Nbi) - r] s (Nbi + 1)} ds\ > 5/2 



Therefore, we just need to prove that, for any 5 > and any t G [0, T] 



lim 

N—f-oc 



J/3 



{r] s (Nh) - ^(iV&i + l)} 



> <5 



= 0. 



Applying Markov's inequality, we bound the previous probability by 



-l E /3 



{?7 S (M> 4 ) - ^(iV&i + l)} 



ds 



Now, in order to conclude it is enough to do the following. First add and sub- 
tract the empirical mean in the box of size eN around Nbi and Nbi + 1- Then, 
by the triangular inequality and since \rfe N (x) — ri E s N (x + 1)| < the term 
involving the two empirical means vanish. For the other two terms, we invoke 
Lemma [5.2[ This finishes the claim. 



^ is concentrated in abso- 



Proposition 6.1. For (3 e [0, 1), any limit point o, 
lutely continuous paths nt(du) = p(t, u) du, with positive density p(t, •) bounded 
by 1, such that p(t, •) is a weak solution of ®. 
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Proof. Let {H l : i > 1} be a countable dense set of functions on C 2 (T), with 
respect to the norm ||i?||oo + ||<9 2 iJj|oo- Provided by d20l ) and intercepting a 
countable number of sets of probability one, is straightforward to extend ( l20l > 
for all functions H G C 2 (T) simultaneously. □ 

6.2. Characterization of Limit Points for /3 = 1. 

The idea in this case is to show that p(t, •) is an integral solution of © for a 
small domain of functions and then extend this set to Hhr. 



Let Cw C % w be the set of functions H in L 2 (T) such that for ieT 



H(x) 



(0,3!] 



for some function h in C(T) and a, 6 G K satisfying 



h{x) dx = 



(0,1] 



Note that a function in is continuous in T\{&i, b k } and well defined ev- 
erywhere. Now, fix a function H g Cw and define the martingale M^(H) as in 
@. We aim that, for every S > 0, the result in ( flit holds for H G Cvv- In fact, 
this was already shown, for H G C 2 (T), in the proof of Proposition ^, ll By ( flOl ). 

for £ G [0,T] 



(M N (H))t<T C+i 



#(#) 



Since if G CV, H is differentiable with bounded derivative, except at the points 
b k . Therefore, for any pair x, x + 1 such that there is no bi between j| 
and ^ii, the following inequality holds 



sc+l ■ 
V ■ 



On the other hand, if there is some {bi : i = 
then £^ x+1 = ^V - ' 3 and in this case we get to: 



1, .., k} in the interval ^jr-), 



h(z&)-h(%) 



< 



N ' ^ N ' 

Since there are only finite k slow bonds, we conclude that the quadratic vari- 
ation of Mf(H) vanishes as N — > oo. Now, Doob's inequality is enough to 
conclude. As above, by a simple change of variables, we may rewrite the mar- 
tingale Mf (H) in terms of the empirical measure as in dl9D . Now we want to 
analyze the integral term in the martingale decomposition d!9t >. 



Lemma 6.2. For any H G C w , 



N- 



0. 
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Proof. Recall the definition of the set Tn given in Section [4] and rewrite the 
previous sum as 



(22) 



Outside bi,...,bk, the operator -^-^ coincides with the Laplacian, and outside 
T N , the discrete operator N 2 h N coincides with the discrete Laplacian. Hence, 
the first term above is equal to 



N E 



N- 



It is easy to verify that H e C 2 (T\{bi, . . . , b k }) and has bounded derivatives. 
Thus, by a Taylor expansion on H, it follows that the previous sum converges 
to zero as TV — s- +oo. On the other hand, the second sum in ( 122D is bounded by 
the sum of 



N E 



dx dW 



xer^ 



and 



E 



N 



Since H e C w , d d 



dx W 



H(^) H(§)\ + N&_ 1>m [H(^) H(§) 

H is a continuous function, therefore bounded. Since Tn 
has k elements, the first sum above converges to zero as N — > +oo. It remains 
to analyze the second sum above, where now the definition of the domain Cw 
is crucial. For each x e T N , one of the conductances above is equal to TV -1 . Let 
us suppose that Cx, x +i = ^ 1 an d £x-i,x = 1> the other case being completely 
analogous. In this case, there exists some bi e ^-]. From the definition of 
Cw and the measure W, the function H has a discontinuity at bi of size 

h(dz) dz . 

Besides that, the function H has also sided-derivatives at 6j of the same value. 
With this in mind, is easy to see that 



H(%)]+N[H(etf)-H(%) 



converges to zero as N — s- oo. Recalling there are finite 2k elements in Tn, we 
finish the proof of the lemma. □ 

Now, fix H € Cw and take a continuous function H £ which coincides with 



H in T\ U** =1 (bi - e, b, + e) and that \\H £ \ 
determined later. Notice that 



< Ili/Hoo. The choice of e will be 



sup \(n t ,H £ - H)\ < sup V/ p(t,u)\H e (u)-H(u)\du<Ake\\H\\ 

0<t<T 0<t<T J (bi-e,bi+e) 
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For every 5 > 0, 



ir. : sup 

0<t<T 



(irt,H) 



> 5 



(23) 



< 



+2Q^ 



7r. : sup 

0<i<T 

7T. : sup 

0<t<T 



^,£^H)d S > 6/3 



{■Kt,H e - H) > 5/3 



By a suitable choice of e, the second probability in the sum above is null. Since 
H e and -£-^H are continuous, by the Portmanteau's Theorem and Proposi- 
tion Ell it holds that 



7r : sup 

0<t<T 



)ds > 5/3 



JV->oo 



lim Pf* 

TV-s-oo 



7r : sup 

0<t<T 



Trt.ff 6 ) - (7r 0! F £ ) - / (ir s ,±^H)ds > 5/3 



sup 

0<t<T 



)ds > 5/3 



Notice that the last equality is just the definition of the measure Since 
there is only one particle per site, it holds that sup 0<t<T \ (tt^,H £ — H)\ < 
4fce||if|| 00 , since H e coincides with H in T\ U* =1 (b l — s,bi + s). Adding and 
subtracting (yrf , N 2 l, N H), (vrf , H) and ,H), we obtain that 



< lim P'L, 



sup \{n?,HZ) (*»,H°) - 

0<t<T 



lim ] 

AT->oo 



-2 lim P£ K 



sup \(ir?,H) - (ttq,H) - / (Trf ,N 2 h N H) ds 

0<t<T I Jo 



> 5/12 



sup (Trf 

L 0<*<T I 



iJ 6 - £T 



> (5/12 



With another suitable choice of e, the third probability in the sum above is null. 
Lemma [6 . 2 1 implie s that the second probability above is zero for N sufficiently 
large. Recall we proved that ( fill holds for H e Cw, so that the first term in 
the sum above is zero. Finally, from the previous computations we conclude 
that ( f23l is zero for any 5 > 0. Therefore, Q* is concentrated on absolutely 
continuous paths n t (du) = p(t, u) du with positive density bounded by 1 and for 
any fixed H G Cw, Q* a.s. 



(p t ,H)-{p ,H) 



Ps,-t^H)ds 



for all f e [0, T] . (24) 
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Proposition 6.3. For [3 = 1, any limit point ofQ^ is concentrated in abso- 
lutely continuous paths nt(du) = p(t,u) du, with positive density p(t 7 •) bounded 
by 1, such that p(t, •) is a weak solution of ©. 

Proof. By a density argument, d24l l also holds, Q* a.s., for all H G Cw simulta- 
neously It remains to extend ( l24l > for H e W-w- For that purpose fix H G Hyy. 
Thus, for xeT 

H(x)=a+ [ (/? + f" ' h(z)dz]w{dy) , 

J(P,x] \ JO ) 

with a,/3 G M, h G L 2 (T) satisfying ®. Let ft„ G C(T) converging to h G L 2 (T). 
Define 



//n(x) = a„ + 



0,a;] 



fin 



h n (z)dz)W(dy), 



where a n -> a and /3„ — > /3. By the Dominated Convergence Theorem, it follows 
that H n converges uniformly to H. Therefore d24t > is true for all H G 'H\ v . □ 

6.3. Characterization of Limit Points for j3 G (1, oo). 

In this regime of the parameter f3, Proposition 15.61 says that is concen- 
trated on trajectories absolutely continuous with respect to the Lebesgue mea- 
sure TT t (du) = p(t,u)du such that, for each interval (b t ,b t+ i), p(t,-) belongs to 
I/ 2 (0,T;'H 1 (6 l ,6j + i)). It is well known that the Sobolev space V}(a,b) has the 
following properties: all its elements are absolutely continuous functions with 
bounded variation, c.f. JH and II 121 . therefore with lateral limits well-defined. 
Such property is inherited by L 2 (0, T; h + i)) in the sense that we can in- 

tegrate in time the lateral limits. Therefore, Qia.s., for each i = l,...,k and 
for any t G [0,T]: 



/ p(s,bf)ds< oo and / p(s,b i+1 )ds 
Jo Jo 



< oo. 



To simplify notation, in this subsection we denote a = and b = b i+ \. Fix 
h G C 2 (T) and define H : [0,T] x T -> M by H(t,u) = h(t,u)l [a>b] (u). 
Recall that ir t (du) = p(t, u)du. We begin by claiming that 

i- 1 



it. : (p t ,H t ) - (p ,H ) 



(p s ,d 2 u H s + d s H s )ds 



d u H{s,a + )p(s,a + )ds+ / d u H(s, b~) p(s, b~) ds = 0, Vi G [0, T] 



= 1. 
(25) 



In order to prove 



it is enough to show that, for every 5 > 



7r : sup 

0<t<T 



(pt,H t ) - (p ,H ) 



(p s ,dlH s +d s H s )ds 



d u H(s,a + ) p(s,a + ) ds + / d u H(s,b ) p(s,b ) ds 



> S 



26 



TERTULIANO FRANCO, PATRICIA GONCALVES, AND ADRIANA NEUMANN 



Since the boundary integrals are not well-defined in the whole Skorohod space 
D([0,T],M), we cannot use directly Portmanteau's Theorem. To avoid this 
technical obstacle, fix e > 0, which will be taken small later. Adding and sub- 
tracting the convolution of p(t, u) with l £ , the probability above is less than or 
equal to the sum of 



7T. : sup 

o<t<r 



(Pt,H t ) - (p ,H ) 



(p s ,dlH s + d s H s )ds 



(26) 



d u H{s,a+)(p s *L e )(a)ds + / d u H(s, b~) (p. * t e )(6 - e) ds > 5/2 



and 



7r : sup 

0<t<T 



d u H(s,a + ) (p s * L £ )(a) ds - / d u H(s,b ) (p s * i e )(b - s) ds 



d u H{s 1 a + )p(s 1 a+)ds+ / d u H{s,b~) p(s,b~) ds > 5/2 



where t e and the convolution p* l £ were defined in d!7l >. The convolutions above 
are suitable averages of p around the boundary points a and b. Therefore, as e 4- 
0, the set inside the previous probability decreases to a set of null probability. 
It remains to deal with ( 126D . 

By Portmanteau's Theorem, Proposition [8]3] and since there is only one par- 
ticle per site, ( 1261 is bounded from above by 



lim 

ri 



7T. : sup 

0<t<T 



n t ,H) - (tt ,H ) - / (TT s ,dlH s + d s H s }ds 



- / d u H(s,a + )(Tr s *L E )(a)ds+ / d u H(s, b~) (tt s * i e )(b - e) ds > 5/2 
Jo Jo 

Now, by the definition of Qf^, we can rewrite the previous expression as 



lim 

N— >oo 
r-t 



sup 

0<t<T 



ir?,dtH s +d s H s )ds 



- / d u H(s,a + )r] e s N (Na + l)ds+ / d u H(s,b~) r/f (Nb) ds > 5/2 



If we consider the discrete torus as embedded in the continuous torus, Na + 1 
is the closest site to the right of a and Nb is the closest site to the left of b. The 
next step is to add and subtract (71^, N 2 h N H) and the previous probability 
becomes now bounded from above by the sum of 



lim 



sup 

0<t<T 



Tif , H t ) - (< , H ) - / (Trf , N h N H s + d s H s ) ds > 5/4 
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and 

lim 

JV-s-oc 



sup 

0<t<T 



n?,N'-L N H a )ds - / (n?,dtH s )ds 



d u H(s,a + )r, £ ^(Na + l)ds + / d u H(s,b~) rf™ (Nb) ds > 5/4 



Repeating similar computations to the ones performed in Section H] we can 
show dllD for a test function H that depends also on time. Therefore the 
first probability above is null. Now we focus on showing that the second 
probability above is null. Recalling the definition of H(s, ■) above, we have 
that H(s, •) is zero outside the interval [a, b]. Besides that, for the set of ver- 
tices {Na + 2, . . . , Nb — 1}, the discrete operator TV 2 h N coincides with the dis- 
crete Laplacian, which applied to H(s, •) converges uniformly to the continuous 
Laplacian of H(s, •). Hence, by the triangular inequality, it is enough to show 
that, for any 6 > 0: 



lim 



sup 

0<*<T 



{N 2 h N H s (^) - d 2 u H s (^)} Vs (Na) ds 



+ ± / {N 2 h N H s (^)-d 2 u H s (^)}r h (Na + l)d S 



i f {N 2 h N H s (^)-d 2 u H s (^)} Vs (Nb)ds 
Jo 

i f {N 2 h N H s (^) d 2 u H s (^)} Vs (Nb+l)ds 
Jo 



d u H(s,a + )rj £N (Na + l)ds + I d u H(s,b~) if s lv (Nb) ds 



> 5 



= 0. 



Since h G C 2 (T), the term involving the Laplacian above is bounded. Now, by 
the triangular inequality, it is sufficient to show that, for any 5 > 0: 



lim I 

N—too 



sup 

0<t<T 



N\, N H s (^)r ls (Na)ds + / NL N H s (*j±±)i ls (Na + l)d 



+ NL N H s (m Vs (Nb)ds+ NL N H s (^±±)r 1s (Nb + l)ds 



d u H(s,a + )r] £lv (Na + l)ds + / d u H( S ,b-)rf s N (Nb) ds 



> 5 



0. 



For each one of the four vertices appearing inside the previous probability, the 
operator Ltv has two conductances, one equals to N^ 13 and the other equals to 1. 
Since /? > 1, the terms involving N~@ converge to zero. The terms involving the 
conductances equal to 1, converge to plus or minus the lateral space derivatives 
of H. Recall from definition of H that d u H(s,a~) = d u H(s,b + ) = for all 
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< s < t. From this, it remains to show that for any S > 



lim 



N- 



sup 

0<t<T 



d u H(s,a + )7 ]s (Na + l)ds- / d u H(s,b~) rj s (Nb) ds 



d u H(s,a + )r) e s ly {Na+l)ds+ / d u H{s,b-) rf s N (Nb) ds 



> S 



is null. Last expression is bounded from above by 



lim 



lim 



J/3 



J/3 



sup 


/ 


L 0<i<T 




sup 




L 0<i<T 





d u H(s,a+){r ls {Na + l)-7 1 £ s N (Na + l)}ds\ > 5/2 
d u H(s,b-){r ls (Nb)-if s N (Nb)}ds\ > S/2 



The integral inside the probability above is a continuous function of the time 
t. Moreover, it has a bounded Lipschitz constant. The same argument as the 
one used in ( 12 11 1 together with Lemma [5.41 imply that the previous expression 
converges to zero when e \, 0, which proves i 



Proposition 6.4. For /3 e (1, oo), any Zimi^ poinZ of{Q^ : N > 1} is concen- 
trated in absolutely continuous paths n t (du) = p(t,u)du, with positive density 
p(t, •) bounded by 1, such that p(t, •) is a weak solution of © in each cylinder 

[0,T] x [bi,b i+1 ]. 



Proof. Given ( 1251 *. it remains to extend the result for all functions H and all 
cylinders [0,T] x [bi,b i+ i] simultaneously. Intercepting a countable number of 
sets of probability one and applying a density argument as in Proposition 16. 11 
the statement follows. □ 



7. Uniqueness of Weak Solutions 

The uniqueness of weak solutions of © is standard and we refer to IflOll 
for a proof. It remains to prove uniqueness of weak solutions of the parabolic 
differential equations © and ©. In both cases, by linearity it suffices to check 
the uniqueness for 7( ) = 0. Notice that existence of weak solutions of ©, 
© and © is guaranteed by tightness of the process as proved in Section [H 
together with the characterization of limit points as proved in Section [6l 

7.1. Uniqueness of weak solutions of ©. 

Let p : R + xT^lbea weak solution of© with 7 = 0. By Definition[3l for 
all H G and all t > 

= [{"-im 11 )*'- (27) 

From Theorem 1 of [5], the operator -^^7 has a countable number of eigen- 
values {A„ : n > 0} and eigenvectors {F n : n > 0}. All eigenvalues have finite 
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multiplicity, = A < Ai < • • • and lim„^oo A„ = oo. Moreover, the eigenvectors 
{F n : n > 0} form a complete orthonormal system in L 2 (T). For t > 0, define 

Notice that R(0) = and since p t belongs to L 2 (T), R(t) is well defined for 
all t > 0. By (|27j>, it follows that f t {p t ,F n ) 2 = -2X n {p t , F n ) 2 . Thus 

dm) — Ej^j-jto,^. 

because n"7i+S„) ^ t ' f ") 2 conver S es uniformly to E„gw n =»(i+A„) f ") 2 » 

as N increases to infinity. Therefore R(t) > and (^i?)(t) < 0, for all t > 
and since R(0) = 0, it follows that R(t) = for all t > 0. As a consequence of 
{F„ : ?i > 0} being a complete orthonormal system, it follows that (p t , pt) = 0, 
which is enough to conclude. 

7.2. Uniqueness of weak solutions of ©. 

At first, we begin with an auxiliary lemma on integration by parts. 

Lemma 7.1. Let p(t,-) be a function in the Sobolev space L' 2 (0, T; V}(a, b)). 
Then, for any H e C^flO, T] x [a, b}): 

p(s, u) d u H(s, u) du ds 

d u p(s,u) H(u,s) duds + / p(s, b) H(s, b) — p(s, a) H(s, a) > ds . 
>o Ja Jo ' 

Notice the partial derivative in p is the weak derivative, while the partial 
derivative in H is the usual one. Besides that, the function H is smooth, but 
possibly not null at the boundary [0, T] x {a, b}, and therefore is not valid the 
integration by parts in the sense of L 2 (0,T;H 1 (a, &)), which has no boundary 
integrals. 

Proof. Fix e > and write H = H e + (H — H e ), where H E coincides with H in 
the region [0, T] x (a + e, b - e), has compact support contained in [0, T] x (a, 6) 
and belongs to C ' 1 ([0, T] x (a, 6)). By the assumptions on H e , we have that 

rT rb 

p(s, u) d u H(s, u) du ds 
= — / d u p(s,u) H s (s,u) duds + / / p(s,u)d u (H — H e )(s,u) duds . 

JO Ja Jo Ja 

Last result is a consequence of H e having compact support strictly contained 
in the open set (a, b). Let f e : [a, b] — >• M be the function such that /(u) = 1 if 
u e (a + e,b — s), /(a) = /(&) = 0, and interpolated linearly otherwise. The 
decomposition H = H f e + H(l - f e ) can be done, but now the function H f s 



J a 
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does not have the properties as required above for H £ . Nevertheless, taking a 
suitable approximating sequence of functions H e , it follows that 

fT P b 

p{s, u) d u H(s, u) du ds 



J a 

T c-b 

J {d u( o(s, u)H(s, u)f{u) + p(s, u)d u (if (s, u)(l - f(u))j } du ds. 
Taking the limit as e I yields the statement of the lemma. □ 

Let p(t, •) be a weak solution of© with 7 = 0. Provided by Lemma [7TT1 for 
any function H e C h2 ([0,T] x (b h b i+1 )), 

Pt(u)H(t,u) du + / / < d u p s (u)d u H(s, u) — p s (u)d s H(s, u) >duds = 0. 
bi Jo Jbi J 

From this point, uniqueness is a particular case of a general result in HUH , 
namely Theorem III. 4.1. In sake of completeness, we sketch an adaptation 
of it to our particular case. Denote by W\ T = VF 2 1 T ([0,T] x (a, b)) the space 
of functions with one weak derivative in space and time, both belonging to 
i 2 ([0,T] x (a, b)) and vanishing at time T. By extending the previous equality 
to H G W\ T it follows that 

J J^ +1 {d uPs (u)d u H(s,u)- p s (u)d s H(s,u)}duds = 0. (28) 
It is not difficult to show that the function 

r 

H(s,u) = — I p(r,u)dr 



belongs to W\ T . Replacing last function in j28l ). then we can rewrite d28l as 

{^d s (d u H(s,u)) 2 -(d s H(s,u)) 2 }duds = 0. 



2 



0. 



10 Jb 

By Fubini's Theorem we get to 

|(cUf(7») 2 - (d u H{0,u)) 2 }du- [ [ ' +1 {d s H{s,u)) 2 duds 

bi L J Jo Jbi 

By the definition of H , its weak space derivative vanishes at time T, so that 
the first integral above is null. Therefore, d s H is identically null, and by the 
definition of H above, this implies that p vanishes, finishing the proof. 

8. Appendix 

Proposition 8.1. Denote by i?Ar(^jv|^a) the entropy of a probability measure 
Pm with respect to a stationary state v a . Then, there exists a finite constant 
K := Kq(o) such that Hjy(pNWa) < K n N , for all probability measures pm- 
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Proof. Recall that v a is Bernoulli product of parameter a. By the explicit for- 
mula given in Theorem Al.8.3 of BIOL 



H N (n N \v a ) = 22, ^Niv) log 



r,e{oa} T iv ay " 
r,e{04} T « 

< Y] VN(V) log r , / yHy" 

' a A (1 — a) \" 

= AT(-log[aA(l-a)]). 



□ 



Proposition 8.2. Assume that L is a reversible generator with respect to an 
invariant measure v in a countable space-state E, and V : M + x E — > E is a 
bounded function. Notice that L + V t will be a symmetric operator in L 2 {v). 
Denote by T t the largest eigenvalue of L + V t : 

r*= sup |<v^ 2 ), + <£/,/>,,} ■ 

Then, the supremum above can be taken over only positive functions f, or else, 

r t = sup {(F t ,(y7) 2 ), + (LV7,A/7),}- 

/ density ^ ^ 

Proof. It follows from the expression of the Dirichlet form (see H1QB ), 

{Lf,f) v = -\ £ v{x)L{x lV )[f{y)-f{x)]\ 

x,y£E 

and the inequality ||/(y)| - |/(ar)|| < \f(y) - f(x)\. □ 

Proposition 8.3. If G\, C7 2 , G% are continuous functions defined in the torus 
T, the application from D([0, T], M.) to K £/ia£ associates to a trajectory {7r f : < 
t < T} the number 



(7T t) Gi) - (7r ,G 2 ) - / (n s ,G 3 )ds 



sup 

0<t<T 

is continuous for the Skorohod metric in D([0, T],M). 

Proof. If G is a continuous function in the torus, the application tt m> (tt, G) 
is a continuous application from .M to R in the weak topology. From this 
observation and the definition of the Skorohod metric as an infimum under 
reparametrizations (c.f. IflOl ). the statement follows. □ 
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